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Adding a non-adsorbing polymer to passive colloids induces an at- 
traction between the particles via the 'depletion' mechanism. High 
enough polymer concentrations lead to phase separation. We com- 
bine experiments, theory and simulations to demonstrate that using 
active colloids (such as motile bacteria) dramatically changes the 
physics of such mixtures. First, significantly stronger inter-particle 
attraction is needed to cause phase separation. Secondly, the fi- 
nite size aggregates formed at lower inter-particle attraction show 
unidirectional rotation. These micro-rotors demonstrate the self as- 
sembly of functional structures using active particles. The angular 
speed of the rotating clusters scales approximately as the inverse 
of their size, which may be understood theoretically by assuming 
that the torques exerted by the outermost bacteria in a cluster add 
up randomly. Our simulations suggest that both the suppression of 
phase separation and the self assembly of rotors are generic features 
of aggregating swimmers, and should therefore occur in a variety of 
biological and synthetic active particle systems. 
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Motile bacteria are simple examples of 'living active mat- 
ter'. Passive particles with sizes similar to those of 
most bacteria, viz., 0.2 — 2/im, are colloids. Such particles are 
in thermal equilibrium with the surrounding solvent, and un- 
dergo Brownian motion. In contrast, self-propelled bacteria 
are active colloids. Such particles function far from equilib- 
rium. This renders their physics far richer than that of passive 
colloids, mainly because they are not subject to thermody- 
namic constraints such as detailed balance or the fluctuation 
dissipation theorem. Thus, bacteria are able to harness their 
activity to power externally-added micro gear wheels [TJ [2] [3] , 
self-concentrate and cluster to form a variety of patterns due 
to geometry, steric effects or biochemical cues ^1 E] • 

Currently, no general statistical mechanical theory relates 
the microscopic properties of individual active particles to the 
macroscopic behavior of large collections of such particles. Re- 
cent experiments on non-interacting suspensions of synthetic 
swimmers [7] show that, as in a dilute suspension of passive 
particles, there is an exponential distribution of particles with 
height, but with an increased sedimentation length. To date, 
however, there has been no experiment designed specifically 
to probe the effect of activity on macroscopic properties that 
arise from inter-particle interaction, such as phase transitions, 
perhaps the quintessential many-body phenomenon. Here, we 
report a systematic study of the physics of phase separation 
and self assembly in a suspension of interacting active col- 
loids in the form of mutually attracting motile bacteria. Our 
experimental results, supported by theory and simulations, 
provide a foundation for general treatments of the statistical 
mechanics of interacting active particles. 

Inter-particle attraction in passive colloids leads to aggre- 
gation and phase separation. Such attraction can be induced 
by non-adsorbing polymers [8]. The exclusion of polymers 
from the space between the surfaces of two nearby particles 
gives rise to a net osmotic pressure pushing them together. 
The range and depth of this 'depletion' attraction is controlled 
by the size and concentration of the polymer, c p , respectively. 



Increasing c p leads to aggregation and, ultimately, phase sep- 
aration. For poly disperse or somewhat non-spherical parti- 
cles, the phase separation is of the 'vapor-liquid' (VL) type, 
with coexisting disordered phases differing in particle concen- 
tration, analogous to vapor and liquid phases in atomic and 
molecular systems. We recently demonstrated that a suspen- 
sion of non-motile Escherichia coli bacteria phase separated in 
this fashion in the presence of non-adsorbing polymers [9] [10] ; 
i.e., non-motile bacteria behave like passive colloids. 

Simple estimates suggest that activity should have a 
strong effect on the depletion-driven aggregation of E. coli. 
The force propelling each cell is ~ 37rr]crv, where a ~ 1/xm 
and v ~ 10 /xm/s are typical cell size and speed respectively 
and 7] rsj 1 cP is the viscosity of the aqueous medium; this gives 
a propulsion force of ~ 0.1 pN. The typical force due to deple- 
tion can be estimated by the ratio of its contact strength at the 
onset of phase separation (~ ksT, the thermal energy) and its 
range (~ the size of the polymer, e.g. as measured by its gyra- 
tion radius, r g ~ 10-100 nm), so that k B T/r g - 0.05-0.5 pN. 
Thus, the competition between activity and inter-particle at- 
traction should generate interesting and novel physics. 

Here we report experiments showing that activity sup- 
presses phase separation: significantly more polymer is needed 
to cause VL coexistence in a suspension of motile E. coli com- 
pared to the passive case [9]^]. A simple calculation suggests 
that this effect can be accounted for quantitatively by an 'ef- 
fective potential' determined by force balance. Simulations 
of active, self-propelled dumbbells subject to a non-specific 
attractive two-body potential support this interpretation. 

Intriguingly, in the range of c p where non-motile cells 
phase separate and motile cells do not, microscopy reveals 
self-propelled and unidirectionally rotating finite clusters of 
cells. Simulations suggest the formation of such 'self assem- 
bled micro-rotors' may be a generic effect in attractive active 
colloids. Our work therefore opens up a novel route to self- 
assembled structures, by exploiting motile activity directly. 

We find that the angular velocities of our self- assembled 
bacterial rotors approximately scale as their inverse size. We 
propose a hydrodynamic theory which explains this finding, 
by assuming that the bacteria on the surface of the clusters 
exert torques that add up randomly to give a unidirectional, 
clockwise or anti-clockwise, rotation. Our simulations provide 
further evidence for this picture, showing that collisions be- 
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A. Cuvettes of E. coli + NaPSS mixtures with the same cell concentration 



(4 X 10 cfu/ml; cfu = colony forming unit) but different amount of polymer (as 
labeled, in wt.%) 2 hours after preparation. Pairs of cuvettes with non-motile (NM) 
and motile (M) cells are shown at each polymer concentration. B. Phase diagram 
showing data from NM (red circles) and M (black squares) cells; filled (open) symbols 
= single phase (phase separated). Approximate positions of the phase boundary are 
indicated by lines. Note that at a cell concentration of 10 10 cfu/ml, the cell bodies 
occupy about 1% of the sample volume. 



tween active particles lead to their random ordering in the 
rotating clusters. 

Note that the only ingredients necessary to explain all our 
observations are self-propulsion and a non-specific short-range 
inter-particle attraction. Thus, the suppression of phase sepa- 
ration and the self-assembly of active rotors should be generic 
features across a wide range of active systems. 



Fig. 2. A. Bright field microscopy image of a motile bacteria-polymer mixture 
(2.4 X 10 9 cells/ml, 0.2% NaPSS). (b) Snapshot of simulation of active dumbbells 
with parameters chosen to be similar to the experiments shown in A. 
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Fig. 3. (A) A schematic of the depletion potentials for a passive (blue) and an 
active (red) system (with the range grossly exaggerated). (B) A typical fitted prob- 
ability distribution function of the swimming speed of our E. coli from DDM [6]. 
The dashed lines indicate the typical velocity values used to estimate F a . (C) The 
calculated phase boundary for a passive (green) and active (red and blue) bacteria- 
polymer mixtures interacting via the depletion potential. The colours correspond to 
the propulsion velocity values in (B) used to estimate the active force F a . 



Results 

Experiments on bacteria-polymer mixtures. We studied sus- 
pensions of motile E. coli bacteria ('smooth swimming' strain 
HCB437 JJ) in motility buffer, and added various concen- 
trations of sodium polystyrene sulfonate (NaPSS, molecular 
weight = 64,700 g/mol). (See Materials and Methods for de- 
tails.) Previous work has shown that NaPSS is non-adsorbing 
to K-12 derived E. coli in motility buffer, and causes depletion- 
driven phase separation at high enough concentrations [9j [10] . 
The range of the depletion attraction is ~ 35 nm [lOj . 

Figure 1A shows samples with identical cell concentrations 
but different NaPSS concentrations (c p ) 2 hours after prepara- 
tion. At each c p we show samples with non-motile (NM) and 
motile (M) cells, where the former were obtained from the 
latter by intense vortexing to break off their flagellaQ Dif- 
ferential dynamic microscopy (DDM) [6] returned an average 
swimming speed of 7 — 10/xm/s for the M cells, and confirmed 
that the NM cells were purely Brownian in their dynamics. 

In the images shown in Fig. 1A, the brightness increases 
with cell concentration. Thus, in both NM and M samples 
with c p = and 0.1%, cells remain homogeneously distributed 
after 2 hours Q While this remains true of the M sample at 
c p = 0.2%, the NM sample in this case has started to separate 
into upper and lower phases with low and high concentrations 
of cells ('vapor' and 'liquid' respectively). This phase separa- 
tion process is complete in the NM sample with c p = 0.5%, 



while it is only just starting at this c p for the M cells. Phase 
separation is nearly complete at the next two higher c p (1% 
and 2%) in the M samples. Similar samples series at other 
cell concentrations enabled us to construct a phase diagram, 
Fig. IB. The phase boundary and phase separation kinetics 
of NM cells correspond to that reported previously [9] ITU] . 
Our first important finding is that the phase boundary for 
the M cells is shifted significantly upwards: more polymer, 
or equivalent ly, a stronger depletion attraction, is needed to 
cause phase separation when the particles are active. 

We next imaged NM and M samples in the 'pre-transition' 
region of the phase diagram, i.e. at polymer concentrations 
just below the respective phase boundaries. In this region 
of a classical atomic or molecular system, we expect to find 
transient clusters (or 'liquid droplets'). We found such pre- 
transition clusters in previous simulations of spherocylinder- 
polymer mixtures designed to mimic our experimental mix- 
tures of NM bacteria and polymer [9] [10] . We observe pre- 



1 Experiments using NM cells of other origin (e.g. killing M cells by heating to 60° C for 1 hour, 
deletion of key flagella synthesis gene FliF) gave the same results. 

2 DDM showed that the motility of cells remained more or less constant over this duration. 
3 Strictly speaking, this is the polymer concentration in a polymer bath in osmotic equilibrium with 
the sample, and not the concentration of polymers in the sample itself; but at the sort of particle 
concentrations we are working at, the two quantities are nearly equal. Under our experimental 
conditions, the polymer concentration around a swimming bacterium is constant. 
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transition clusters in both NM and M samples: the static 
snapshots are rather similar (an example from an M sample 
is shown in Fig. 2a). Movies of the clusters, however, reveal a 
dramatic difference between NM and M clusters. NM clusters 
undergo what is recognizably Brownian motion in both trans- 
lational and rotation (Supplementary movie SI). Remarkably, 
however, the clusters self- assembled from motile cells translate 
with speeds approaching that of single cells and also rotate 
(Supplementary movies S2 and S3). The clusters persist over 
the entire duration of observation (minutes) , during which the 
rotation is essentially unidirectional in a frame of reference in 
which the axis of the rotation is stationary (see Supplemen- 
tary movies S2 and S3), although the axis of rotation itself 
seems to drift randomly. Occasional interruptions in this uni- 
directional rotation do occur, apparently due to collision with 
single motile cells or other clusters. The sense of the rotation 
appears to be random. This, then, is our second important 
finding. The interplay of motility and depletion attraction 
leads to the self-assembly of spontaneously rotating clusters, 
or micro-rotors. 

Note that our self-assembled micro-rotors are qualita- 
tively different from the bacteria-driven rotation of exter- 
nally added micro-lithographed gear wheels with a built-in 
chirality [TJ [2 [3]. Our micro-rotors are not only powered 
by bacteria, but self-assembled by and from them. The only 
externally-added agents are non-adsorbing polymers. 

Effective potential theory for active phase separation. We 

first propose a phenomenological theory for the observed sup- 
pression of phase separation. Our theoretical framework 
should be applicable to phase separation caused by inter- 
particle attraction of any origin in a system of self-propelled 
particles of any shape. For analytical tractability, we model 
individual bacteria as hard spheres of diameter a in the pres- 
ence of non-adsorbing polymers. The latter are modeled as 
mutually inter-penetrable spheres of diameter £cr that are 
impenetrable relative to the hard spheres. In this Asakura- 
Oosawa (AO) model for colloid-polymer mixtures, the deple- 
tion attraction between two spheres at £ <C 1 can be approx- 
imated by [4]: 

{oo r < a, 

a<r<a(l + 0, [1] 
r><7(l + 0- 

Here e — ^4> p kBT (1 + / £ is the attractive well depth, and 
4> p is the volume fraction of polymer coils [^] 

In the absence of activity, two particles stay bound at sep- 
aration a until a thermal fluctuation increases this distance 
to beyond a (1 + £), at which point the particles cease to in- 
teract and become free. To estimate the effect of activity on 
depletion-driven phase separation, we consider two active par- 
ticles pushing in the opposite directions with an active force 
±F a while being held stationary by the depletion force. The 
particles come free at separation r* , where the total force act- 
ing on the particles becomes zero: 

-^ + F„ = 0, [„ 

giving r*/a = 1 + f - £k b Tf/2e, with / = F a ^a/k B T. We 
take the effective potential confining two active particles to be 

UeB(r) = U(r)-F a r + Uo, [3] 

where Uo is chosen such that U e R (r* ) =0 (see Supplementary 
Material), Fig. 3A. Compared to the bare depletion potential, 
U e s(r) is shallower, and has slightly shorter range. 

To estimate the phase boundary for particles interacting 
via £/eff(r), we use the law of 'extended corresponding states' 



[T], which states that the vapor- liquid phase boundaries (or 
'binodals') of various systems of attractive particles collapse 
onto a single master curve if the attraction is characterized 
by its reduced second virial coefficient, 62 = B2/B^ s , where 
B2 and Bf s are the second virial coefficients of the attractive 
particles and of equivalent hard spheres respectively. The law 
is known to work, e.g., in relating the phase behavior of pro- 
tein solutions and colloid-polymer mixtures [15] . We therefore 
estimate the binodal in our system by mapping onto that of 
Baxter 'adhesive hard spheres' [3] — hard spheres with an in- 
finitely short range attraction (see Supplementary Material). 

We find that the shift in the VL binodal strongly depends 
on the value of the active force F a with which a stationary bac- 
terium is pushing against the depletion potential. This force 
is approximately twice the force produced by the flagella of a 
free-swimming bacterium (see [22] and Supplementary Mate- 
rial for more details), which is proportional to the swimming 
speed, v, of the bacteria. The measured probability density 
of the latter, P(v), can be fitted by a Schulz distribution [6], 
see Fig. 3B, where we indicate by dashed lines the range of 
speeds used for estimating the value of F a . 

In Fig. 3C we plot the contact strength of the depletion 
potential e/ksT required to trigger phase separation for the 
passive and active bacteria-polymer mixtures as a function of 
cell volume fraction. We predict that activity shifts the bin- 
odal upwards by a factor of ~ 1.4—1.9, depending on the 
swimming speed used. Since e/ksT is proportional to c p in 
our system, this should translate into a corresponding shift in 
the experimental phase boundary. 

The predicted shift of the phase boundary by a factor of 
< 2 upwards has the same sign and a similar value as the 
observed upwards shift by a factor of < 3. The latter is likely 
an overestimate, since sedimentation of pre- transit ion clus- 
ters renders the experimental phase boundary for NM cells a 
lower bound. (The motility of the active pre-transition clus- 
ters means that sedimentation has minimal effect on the ob- 
served phase boundary for M cells.) Considering the crudeness 
of the approximations involved in our theory, we take such a 
degree of correspondence between experiment and theory as 
a validation of the basic physical content of our approach, 
viz., that activity renders it easier for thermal fluctuations to 
'free' cells bound by interparticle attraction, an effect that we 
model by a shallower effective potential, Fig. 3A. 

Note that we have above neglected hydrodynamic forces. 
These create an attraction between two particles as they move 
apart, which is proportional to their relative velocity and 
therefore vanishes at the balance point between the active 
and polymer-induced forces. Moreover, hydrodynamic inter- 
actions have no effect on barrier heights or phase boundaries 
in equilibrium systems. Their neglect therefore appears justi- 
fied within our quasi-equilibrium evalation of the interaction 
potential. 



Simulating active phase separation. The simple physical in- 
gredients entering our effective potential theory suggest that 
the suppression of phase separation by activity should be 
generic. To test this, we performed Brownian dynamics (BD) 
simulations intentionally preserving only what we think are 
the essential features of our experiments: motility and a short- 
range attraction (depth e) but without hydrodynamics. Specif- 
ically, we simulated a collection of Brownian dumbbells (each 
with local frictional drag in the lab frame) interacting via a 
non-specific short-range attraction; a constant force applied 
along the axis of each dumbbell renders it motile (see Materi- 
als and Methods and Supplementary Material) . While numer- 
ical parameters were chosen to resemble our experiments, our 
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Fig. 4. The energy of a gel in a simulated system of attractive dumbbells measured 
in units of the contact potential, U/e, as the temperature, measured as e/kgT, 
is changed, for active (bold) and non-active (thin) particles. The energy abruptly 
increases as the gel melts. 
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Fig. 5. Plot of the angular velocity as a function of cluster radius. Different 
kinds of points refer to independent experimental data sets. The green, red and blue 
solid lines are predictions of Eqs. 3 and 4 with F a being 0.05, 0.1 and 0.15pN, 
correspondingly. 



qualitative conclusions are valid for a much wider parameter 
space. 

Our simulated system is too small to allow reliable direct 
estimate of the effect of activity on the position of the binodal. 
Instead, we study how activity renormalizes the e/ksT axis 
by quenching otherwise identical starting systems of active 
and passive attractive dumbbells to low temperatures (large 
e/ksT), at densities where the system forms a space-spanning 
gel. Starting from this configuration, we performed runs at 
decreasing e to locate gel melting as the point at which the 
potential energy per particle, U, abruptly increased, Fig. 4. 
Intriguingly, the melting transition in the active system ap- 
pears sharper than in a passive one. More importantly for our 
purposes, motility increases the value of e which is needed to 
stabilize the gel by a factor of ~ 3. Choosing different initial 
conditions and densities does not much affect this result (see 
Supplementary Materials for more details). 

The substantial agreement between experiment, Fig. IB, 
theory, Fig. 3C, and simulations, Fig. 4, on the effect of 
activity in renormalizing the e/ksT axis, and the absence of 
system-specific features in either our theory or simulations, 
suggest that phase separation suppression by activity is likely 
generic for attractive particles. 



Self-assembled rotors. We have observed coherently-rotating 
pre-transition active clusters in both experiments (Supple- 
mentary movies S2 and S3) and simulations (Supplementary 
movie S4). The existence of pre-transition clusters per se is 
unsurprising. However, their persistence (over at least min- 
utes in experiments) and unidirectional rotation require ex- 
planation. 

First, clusters of short-range attractive particles are long- 
lived due to entropic reasons (Willem Kegel, private commu- 
nication). Clustering in atomic systems interacting via van 
der Waals attraction is expected at a well depth of e ~ ksT. 
A shorter range attraction entails higher loss of entropy upon 
bonding, which therefore requires a larger e. Cluster life time 
is controlled by the escape rate of a single cell bonded to 
its neighbors, which itself is dominated by a Boltzmann fac- 
tor, e e l kuT . A quantitative estimate using the Kramers for- 
malism (see Supplementary Material) returns a life time of 
many minutes to months (depending on the average number 
of neighbors per cell we use in the calculation) , consistent the 
observation of persistent clusters. 

Secondly, since we found self- assembled rotors in our BD 
simulations, which did not include fluid-mediated interactions, 
we may rule out specifically hydrodynamic explanations for 



the phenomenon, and the self assembly of micro-rotors is likely 
generic. We now show that the quantitative features of the 
micro-rotor motion do depend on system-specific details. 

The angular speed, of our rotating bacterial clusters 
(see Materials and Methods) decreases with cluster size, R, 
taken as half the arithmetic average of the longest and short- 
est dimensions of each cluster, Fig. 5. The data is consistent 
with Q rsj BT 1 5 with typical Q in the range 1-20 rad/s. (See 
Supplementary Material for a discussion of the range of cluster 
sizes represented in our data.) 

The physical origin of cluster rotation is clear. The forces 
generated by the flagella bundle of each bacterium (which 
hereafter we refer to as a single (effective) flagellum) in a 
cluster do not cancel, giving a residual net torque about the 
cluster center. We assume that only bacteria on the surface 
of the cluster exert torques, implying that the total number 
of bacteria, N, participating in propulsion and rotation of 
the cluster is equal to N = 4ttR 2 /Ao, with each bacterium 
occupying Aq ~ 2 /xm 2 . If each bacterium exerts a torque of 
magnitude To on the cluster and we assume that these torques 
are randomly oriented, then the total torque on the cluster, 
Ttotj is given by T to t ~ y^NTo, since the sum of N random 
vectors of equal length scales as ~ R. T to t is balanced by 
the rotational friction, which we take to be 87r?7i? 3 Q, where 
rj ~ 1 cP is the solvent viscosity. Thus, we predict that 



:T . 



[4] 



If we take each bacterium to be a force monopole of magni- 
tude F a that contributes a torque To ~ F a R, then Eq. [4] gives 
Q ~ R~ x , which is close to our observations. 

However, the total force exerted by a swimming bacterium 
on the fluid must sum up to zero, so that the lowest or- 
der approximation with the appropriate symmetry is not a 
monopole, but a dipole [171 118j . For simplicity, we assume 
that bacterial flagella lie tangentially to the cluster surface. A 
small portion of the flagellum at some distance from the clus- 
ter generates a propulsion force that is transmitted through 
the whole flagellum to the cluster, while a force of the same 
magnitude and opposite direction is applied to the fluid lo- 
cally. The former corresponds to the monopole contribution 
already discussed. The latter force generates a fluid flow that 
has to vanish at the surface of the cluster, which exerts an ex- 
tra drag force acting on the cluster that partially cancels the 
'direct' force transmitted through the flagellum. Summing 
up the contributions from different parts of the flagellum we 
find (see Supplementary Material) that the torque exerted per 
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bacterium on the surface is 

To = F a R [1 - R(l 2 + R 2 )- - 5 ] , [5] 

with / ~ 10/im being the length of a flagellum. The right hand 
side consists of a monopole term F a R and a dipolar correc- 
tion. We estimate F a as the force necessary to propel a 1/mi 
diameter sphere at speeds in the range 5-15/xm/s, Fig. 3B, 
giving F a ~ 0.05 - 0.15 pN. The prediction of Eqs. [3] and [5] 
for various values of F a , Fig. 5, is compatible with the data. 
Note that when R < /, as in our experiments (Fig. 5), the 
dipolar correction in Eq. [5] does not effectively change the 
scaling, so that we may expect Q ~ R -1 . 

We have assumed that the forces in a cluster add randomly 
and that their relative disposition is fixed. The alternative as- 
sumption that they add coherently predicts an angular veloc- 
ity independent of cluster size, clearly at odds with our data 
in Fig. 5. While it is non-trivial to test experimentally our 
'quenched disorder' assumption, it is supported by an analysis 
of the ordering of active dumbbells within a rotating cluster 
in our BD simulations. These show that orientational order- 
ing is largely absent, and that there is little rearrangement 
of the dumbbells after the clusters form (see Supplementary 
Material) . The neglect of hydrodynamics in our simulations is 
not a major shortcoming in this respect, as in the micro-rotor 
phase we expect intra-cluster interactions to be dominated by 
excluded volume, rather than hydrodynamics. 

Discussion and Conclusions 

We have studied phase transition and pre-transition clustering 
in a mixture of bacteria and non-adsorbing polymers, viewed 
as an active suspensions with inter-particle attraction. We 
find that activity suppresses phase separation compared to 
the case of non-motile bacteria [9 , 10 . Addition of ^ 3 times 
more polymer was needed to cause phase separation in a sus- 
pension of motile cells. We rationalize this by modeling the 
interaction between self-propelled cells by a shallower effective 
potential, determined by balancing the depletion force and the 
self propulsion. Quantitative prediction of this upward shift 
of the phase boundary follows from the extended law of cor- 
responding states [1] via mapping second virial coefficients. 

An alternative interpretation of our results is to ascribe a 
higher effective temperature, T e ff « 3T, to our active particle 
system while leaving the inter-particle potential unchanged. 
The idea of a higher effective temperature due to activity 
has recently been applied to characterise the transport prop- 
erties of a suspension of synthetic self-propelled colloids [7]. 
Each active particle undergoes a random walk in the long- 
time limit because the direction of its self-propelled motion 
is subject to thermal fluctuations (controlled by the thermo- 
dynamic temperature T) [19]. This random walk is described 
by an effective diffusion coefficient D e ff = v 2 tr/4, where v is 
the self-propulsion speed and tr = 4L 2 /3Do is the (Brown- 
ian) rotational relaxation time of an object of radius L with 
(Brownian) translational diffusion coefficient Do. Palacci et 
al. found that D e ff = ksT e ^ /67rr]L. For our smooth swim- 
ming E. coli strain, we may estimate D e ff < 10 4 /mi 2 /s [20] {^] 
giving T eff - 10 4 T. ' ' 

This latter effective temperature for characterizing diffu- 
sive transport in the bacterial suspension is four orders of 
magnitude higher than the effective temperature that we may 
use to characterize a suspension of the same motile bacteria 
in the context of phase separation driven by inter-particle at- 
traction. This 'discrepancy' highlights an important point: 
there is no single parameter, 'the effective temperature', that 
is appropriate for systems away from equilibrium. Indeed, the 
whole concept of effective temperature may be more appropri- 



ate for some situations than others. In a system such as ours, 
we suggest that it is more appropriate to talk of an effective 
potential between particles, which is explored by fluctuations 
that are characterized by the thermodynamic temperature. 

Remarkably, when the inter-particle attraction is just too 
weak to cause phase separation, we observe finite clusters that 
individually break chiral symmetry and rotate unidirection- 
ally. The spontaneous formation of these micro-rotors consti- 
tutes a clear demonstration of the self assembly of functional 
structures from active suspensions. We accounted for the ob- 
served scaling of cluster rotation speed with cluster size using 
a hydrodynamic theory. These 'active clusters' show minimal 
sedimentation: their motion confers a high effective tempera- 
ture in the diffusive sense [7]. 

We have also simulated Brownian dumbbells each subject 
to a body force directed along its axis and with a short-range 
inter-particle attraction, with parameters chosen to mimic 
our experimental bacteria-polymer mixtures. Our theoretical 
framework is also able to account for the melting of gel states 
in the simulated system, which we take to be a surrogate for 
the phase boundary in our small simulations. Moreover, clus- 
ters in the melted state also show spontaneous rotation. 

Our experimental system is a very specific one: living 
bacteria propelled by rotating flagella bundles operating via 
low- Reynolds number hydrodynamics in a chemically-complex 
buffer with added non-adsorbing polymer; the latter causes a 
short-range inter-particle attraction via the depletion mech- 
anism. Our phenomenological theory and simulations strip 
away many of these system-specific features. The success of 
the theory and simulations in reproducing experimentally ob- 
served phase separation suppression and the self assembly of 
rotating clusters suggest strongly that these are generic fea- 
tures in attractive active particle systems. It would therefore 
be interesting to look for such effects not only in suspensions of 
other bacteria, but also in suspensions of self-propelled 'Janus' 
particles [19]. In the latter case, the competition for 'fuel' 
among neighboring particles may act to limit cluster sizes. 

Finally, we point out that our bacterial suspensions con- 
tain a fraction of non-motile organisms [6 . In the presence of 
polymer, these cells should be the first to aggregate, producing 
'nuclei' for the subsequent aggregation of motile cells. Prelim- 
inary BD simulations of mixtures of motile and non-motile at- 
tractive dumbbells support this picture, which motivates our 
assumption that cluster rotation and translation are powered 
by motile cells on the outside. It also points to possible rel- 
evance of our work for aquatic and marine ecology: bacteria 
aggregating in the presence of non-living organic particulate 
matter (detritus) can trap the latter and keep it suspended 
for the water column food chain |23| . 



Materials and Methods 

Culture conditions for motile bacteria. We used E. coli strain 
HCB437, a smooth-swimming mutant. Overnight cultures 
(~ 18 h) were grown in 10 ml Luria Broth (LB, tryptone 
10.0 g/1, yeast extract 5.0 g/1, NaCl 5.0 g/1) starting from 
a single colony on LB agar (tryptone 10.0 g/1, yeast extract 
5.0 g/1, NaCl 5.0 g/1, agar 15 g/1) using an orbital shaker at 
30° C and 200 rpm. A fresh culture was inoculated as 1:100 
dilution of overnight grown cells in 35ml tryptone broth (TB, 
tryptone 10.0g/l, NaCl 5.0g/l) and grown for 4 h (late expo- 
nential phase) . These cells were treated with the following two 
protocols to achieve motile and non-motile cells, respectively. 



4 Note that we have taken 2L ~ 10/j,m to be an estimate of the 'cell body + flagella bundle' 
unit; this is the size that controls Brownian rotation. 
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Preparation of motile cells. Cells were washed three times 
with motility buffer (MB, pH = 7.0) containing 6.2 mM 
K2HPO4, 3.8 mM KH2PO4, 67 mM NaCl and 0.1 mM EDTA 
by careful filtration (0.45 fim HATF filter, Millipore) to min- 
imize flagellar damage. The final volume of the washed bac- 
terial samples was 1-2 ml therefore suspensions could be con- 
centrated up to 25-fold compared to the original culture. 

Preparation of non-motile cells. Cells were washed by cen- 
trifugation (10 min, 2700g Hermle Z323K). After removing 
the supernatant cell pellets were re-suspended in 1 ml MB by 
vortexing for 2 min. In total, three washing steps were com- 
pleted, with the final sample being suspended in 1 ml MB. 

Characterization of bacterial samples. Optical density (OD) 
measurements at 600 nm (Cary IE, Varian) normalized by 
viable plate counts on LB agar of serial diluted samples 
(1.55 x 10 9 cfu/ml OD600nm 1) were used to determine 
cell densities. Motility was characterized by DDM [6]. 

Polymer. Sodium polystyrene sulfonate (Aldrich) was used 
(NaPSS, M w = 64,700 g/mol, M w /M n = 3.1). The molec- 
ular weights and polydispersities were determined by gel- 
permeation chromatography (GPC) against PSS standards. 
We previously estimated a radius of gyration in MB of 17.5 nm 
[9j [TO] . Polymer stock solutions were prepared at 20% (w/v) 
in MB and filtered through a 0.2 /am disposable syringe filter 
prior to use. 

Study of phase separation. Observations were made in closed 
1.6 ml disposable cuvettes with a total sample volume of 1 ml. 
Polymer solution, cell suspensions and MB were mixed in dif- 
ferent ratios to achieve cell concentrations in the range of 10 9 
to 10 11 cfu/ml and polymer concentrations in the range of to 
2 wt%. Samples were homogenized by thorough, careful mix- 
ing, placed inside an incubator at 22C (MIR- 153, Sanyo) and 
observed using a camera (Qlmaging, Micro-publisher 3.3RTV) 
controlled by QCapture pro 5.0 software. Images were cap- 
tured for varying periods up to 24 h. 



Rotation of motile clusters. Two types of sample cells were 
used: capillary cells (dimensions 8 x 50 x 0.4mm, CMS) com- 
pletely filled with 190/il sample (sealed with Vaseline), and 
8-well chambered microscopic cover glass cells (dimensions 
8 x 8 x 8mm, Lab-Tek, Nunc) half-filled with 400/il test sam- 
ple covered with a lid. Samples were prepared at a cell con- 
centration OD600 =0.3 (5 x 10 8 cells/ml) and lwt% NaPSS. 
Clustering of bacteria was observed using either bright-field 
or phase- contrast microscopy at 90 x magnification using two 
different microscope-camera combinations: 1) Nikon Eclipse 
TE2000-U inverted microscope and Marlin F145B2 camera 
(17 fps) and 2) Nikon Eclipse Ti inverted microscope and 
Cool-Snap HQ2 camera (11 fps). Videos were recorded for 
30 to 180 s. Clusters were observed during a period when the 
axis of rotation was perpendicular to the imaging plane. The 
angle of rotation was measured using ImageJ angle tool. Mea- 
surements were made 50-150/xm above the surface at different 
times during the first 3 h after polymer addition. 

Dumbbell simulations. We simulated a system of 1000 dumb- 
bells, each composed by two spheres of radius a and kept to- 
gether by a stiff harmonic spring. Simulations were performed 
by using the LAMMPS code 24 in the (overdamped) Brow- 
nian dynamics mode, and we employed the following addi- 
tional force fields. Depletion-induced attraction was modelled 
by means of a truncated and shifted Lennard- Jones potential, 
where the depth and the size are matched to the interaction 
strength and range of NaPSS-induced depletion. Motility was 
achieved by applying a force directed from the rear to the 
front bead of each of the dumbbell, and applied in its centre 
of mass. Finally, the random force in the Brownian dynamics 
was chosen so as to satisfy the fluctuation-dissipation theo- 
rem in the passive limit. More details of the methods and 
parameter values are given in the Supplementary Material. 
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Phase separation of active colloids 

Here we explain how we estimate the phase boundary for ac- 
tive particles interacting via the depletion potential. First, we 
incorporate the effect of activity by defining an 'effective po- 
tential' based on force balance. Then we invoke the extended 
law of corresponding states [T] and match the second virial 
coefficients along the phase boundary of the parent depletion 
potential and of our effective potential via the known phase 
boundary of Baxter's adhesive hard spheres [2j[3]. The result 
will show that activity significantly postpones the transition. 

Depletion potential between active particles. We model bac- 
teria as hard spheres of diameter a. The interaction potential 
between two hard spheres in the presence of non-adsorbing 
polymer of diameter £cr with £ <C 1 can be approximated by 



U(r 



Here, 



r < a, 

[j-l-£] 2 *<r<a(l + 0, [1] 

r ><r (! + £)■ 



[2] 



is the magnitude of the potential at contact, U(a) — — e, 
and r/p is the polymer volume fraction^] In the absence fluc- 
tuations, the lowest-energy configuration for two particles is 
touching, i.e. at a center-to-center separation of a. 

We model an active particle as a sphere subject to 
a propulsive force, and quantify the competition between 
polymer-induced attraction and activity for the case of two 
active particles with equal and opposite propulsive forces 
The force generated by the depletion potential is given by 



F(r) 



dU{r) 
dr 



2e fr 



i 



[3] 



and is maximal at contact: F c — — 2e/£cr. To cause phase 
separation, e > ksT. Our polymer diameter is £cr ~ 35nm, so 
that F c > 0.2 pN. The propulsive force generated by a free- 
swimming bacterium can be estimated by the drag force on the 
cell body, Fprop °° ^rjav with rj = 10 -3 Pa-s, v ~ 10/xm/s, 
and a ~ 1/mi, so that Fprop ~ O.lpN. These have similar 
orders of magnitude, so that we expect activity to be a signif- 
icant perturbation on phase separation. 

Two passive particles stay bound at distance a until a 
thermal fluctuation increases the inter-particle separation be- 
yond a (1 + £), at which point the particles become free. To 
see how activity changes this scenario, we first note that the 
total force on each active particle is given by the sum of the 
depletion and active forces: 



F eff (r) =F(r) + F p5 



[4] 



This bound pair breaks up when a thermal fluctuation in- 
creases the inter-particle distance to some value r* such that 
i ? e ff(r*) = 0, which, from Eq. [2], is 

^k B T 
2" 



where / = F prop ^a /ksT is the typical work done by the ac- 
tive force in separating the particles by the range of the deple- 
tion potential in units of the thermal energy. Note that since 
/ oc Fprop and Fprop oc v, f increases linearly with swimming 
speed. 



<J 1 


r* (7(1 + 6 











Fig. SI. Potential of interaction between colloidal particles. Blue: depletion po- 
tential U(r) in the absence of activity; red: sketch of the effective potential U e ff(r) 
for self-motile particles; r* marks the point where the total binding force due to de- 
pletion and activity vanishes. 



We obtain the effective potential of interaction, L r e ff(^), 
between two particles in the presence of depletion and activ- 
ity by requiring that 



dr 



which after integration gives 

Ues(r) = U(r)-F prop r + U , 



[6] 



[7] 



where the integration constant Uo is chosen so that U e & (r* ) = 
0, since the particles become free after being separated beyond 
the distance r* . This procedure yields 







r — r* 



r < a, 

a < r < r*, 

r > r*. 



In Fig. SI we sketch the bare depletion potential U(r) in the 
absence of activity and the 'active' depletion potential U e ft(r). 
We observe that the active force decreases both the range of 
the potential and its slope. Note that this procedure is widely 
applied in modeling the effect of an applied force on ligand- 
receptor binding [5]. 



Reserved for Publication Footnotes 



1 Strictly, this is the polymer volume fraction in a reservoir of pure polymers in osmotic equilib- 
rium with the bacteria-polymer mixture. However, at the bacterial concentrations we work at, this 
distinction is unimportant. 

2 Note that any other orientation of active forces is unstable and will convert into this one, with the 
active forces pointing in opposite direction. 
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The second virial coefficient for the 'active' depletion po- 
tential is given by 



B. 



(adp) 
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Here 
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where we have introduced x = e/ksT. 



e dt, [11] 



Baxter potential. The Baxter adhesive hard sphere potential 
[2 is defined by 



U B (r) 
k B T 



^ [I2r(^)] 




r < a, 



r > , 



[121 



in the limit d — »• a. Here d — a is the range of the attrac- 
tive potential and the Baxter effective temperature r sets its 
strength. Its second virial coefficient is: 



B. 
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[131 



The phase diagram for this system was calculated by Miller 
and Frenkel [3]. 

The extended law of corresponding states. We map a system 
of active particles interacting via the depletion potential onto 
the Baxter model by matching their second virial coefficients 
along the gas-liquid binodal. This mapping yields a relation 
between the Baxter effective temperature r and the parame- 
ters of the active suspension. 



Solving B^ p) 



7-,(adp) 
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[141 



This expression is used to calculate the phase boundary for 
the active depletion potential. For any value of the particle 
volume fraction 0, we use the position of the phase boundary, 
r ((/>), found by Miller and Frenkel 3 and find from Eq. |l4| 
the corresponding depth of the depletion potential x = e/ksT 
at phase separation. We repeat this procedure for different 
values of the dimensionless active force / to study how activ- 
ity influences the phase diagram. 

Estimating the propulsion force. Here we provide a more de- 
tailed estimate for the propulsion force Fprop. Note first that 
there is a distribution of swimming speeds in our bacterial 
populations. A representative swimming speed distribution 
from fitted differential dynamic microscopy [6 data is shown 
in Fig. 3B in the main text. Since the average speed in this 
population is about v ~ 10/mi/s, a significant fraction of or- 
ganisms swim at > 10/xm/s. 

Within the resistance matrix framework of a free swim- 
ming E. coli [3|8], the propulsion force Fprop generated by 
the flagella bundle and the motor torque N m required to ro- 
tate it are expressed in terms of the propulsion velocity v and 
the rotation speed of the flagella bundle uj: 

Fprop = -Av + Buj, [15] 
N m = -Bv + Du, [161 




Fig. S2. The depth of the depletion potential at which phase separation occurs as 
a function of particle volume fraction, for different values of activity (from bottom to 
top): green: passive suspension (/ = 0), red: / = 2 (corresponding to the average 
speed of V = 10/im/s), blue: / = 5 (corresponding to v = 20/im/s, close to the 
maximum speed observed in our bacterial populations). The shaded area therefore 
represents the shifts that can be expected from the faster half of our cell population. 



where A, B and D are the resistance matrix coefficients re- 
lated to the bundle geometry. For a non-tumbling strain of 
E. coli swimming with v — 20/xm/s and uj — 780rad/s, Chat- 
topadhyay et al. [9] found that 



A = 1.48 x 10" 8 Ns/m 
B = 7.9 x 10" 16 Ns, 
D = 7.0 x 10" 22 Nsm 



[17] 



For a free swimmer, these values give Fprop ~ 0.3pN, which is 
slightly higher than the estimate used above. 

Now consider a bacterium with free-swimming speed v — 
20/xm/s that is held stationary by a depletion potential. It is 
easily verified that for all v < 20/xm/s, Bv< Duo in Eq. ( 16 ). 
Thus, for our stationary bacterium with v = 0, iV m , and hence 
uj, remain approximately equal to their free-swimming values. 
On the other hand, the propulsion force, Eq. (15), is signif- 
icantly increased to Fprop = B uj ~ 0.6pN, leading to / ~ 5. 
In what follows, we study the effect of activity on phase sep- 
aration for / < 5. 



Phase diagram for the active depletion potential. Our results 
are plotted in Fig. S2, which corresponds to Fig. 3C in the 
main text. It is clear that activity significantly suppresses 
phase separation, by up to a factor of ~ 1.9, which corre- 
sponds to / = 5 (or a swimming speed of v = 20/im/s). 

Note that we have used a similar approach previously 
to explain how oscillatory shear drives the crystallization of 
depletion-induced colloidal gels [10] . Consider two particles 
bound by depletion in an initial configuration such that at 
the extremes of each oscillatory cycle, these two particles are 
maximally separated by the imposed shear. In this new con- 
figuration, the particles experience a shallower effective po- 
tential. Estimating the particles Kramers escape time out of 
this shallower effective potential provided a semi-quantitative 
explanation for the frequency and amplitude dependence of 
shear-induced crystallization. The details of the treatment 
differ because [10] involves an imposed strain, and the present 
work involves an imposed stress. 
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Cluster life time 

Here we expand on our explanation for the persistence of pre- 
transition active clusters. A Brownian particle (diffusion co- 
efficient D) trapped in a depletion potential of depth e and 
range 5 approximated as a ramp of the same depth and range 
has a Kramers escape time given by [10] 



t K 



6Dx 2 



[181 



where x = e/ksT. Differential dynamic microscopy [6 mea- 
sures D & 0.3//m 2 /s for our bacterium. The range is esti- 
mated by 5 — r* — a, Fig. SI and Eq. |5|, which, like x, is a 
function of the propulsion speed via the oimensionless param- 
eter /. At the average swimming speed of v = 10/xm/s, / = 2, 
from which we obtain 8 ~ 30nm (using a — 35nm, which is 
the size of our polymer) and x ~ 7 at the phase boundary. 
Each cell on the surface of the cluster is in contact with n 
other cells, so that its escape time, and therefore a 'cluster 
life time' can be estimate using Eq. ( 18 ) with x replaced by 
nx. Using n — 3 we obtain Ik ~ 1500s (or ~ 25 minutes); 
if we take n — 4, this increases to tx ~ 26 x 10 6 s (or just 
under 10 months). These long escape times explain why our 
pre-transition clusters do not appear to break up and reform. 



Brownian Dynamics (BD) simulations 

Here we detail the model and methods used to simulate a sus- 
pension of active, self-propelled, Brownian dumbbells, discuss 
the mapping from simulation parameters to experiments, and 
spell out some technical aspects of the data analysis. 

Simulation details. In our BD simulations, we study the evolu- 
tion in three-dimensional space of a system of N = 1000 hard- 
dumbbells in an NVT ensemble using cubic periodic boundary 
conditions. A dumbbell consists of two spheres with the same 
diameter a. Any pair of spheres in the simulation interact via 
a truncated and shifted Lennard- Jones potential 



V LJ (r) 



4e 



(?)"-(;)' 



[191 



for r < r c [r — n — rj, with i,j — l,...,2iV), whereas 
Vlj{t) — for r > r c . Choosing a cutoff of r c = 1.2cr gives a 
very short range attraction that mimicks the depletion poten- 
tial. The strength of Vlj is controlled by e, which is varied 
to simulate a change in the concentration of the polymer in- 
ducing the depletion attraction. Finally, the two spheres in 
a dumbbell are 'glued' together by means of a stiff harmonic 
potential Vh(t) — n{r — a) 2 , where k — 700/csT/cr 2 is the 
spring constant. 

Our simulations are carried out using the open source 
LAMMPS Molecular Dynamics package [11 . The motion 
of the spheres are governed by the following under-damped 
Langevin equations of motion: 



dt 2 



-c 



dr* _ dV_ 
dt dvi 



+ F r + F a 



[201 



where m is the mass of a sphere, £ is the friction coefficient 
(C = m7 with damping coefficient 7), V is the total conserva- 
tive potential acting on each particle (V = Vlj + Vh) and F r 
the random force due to the solvent at temperature T. Ac- 
tivity via self-propulsion is introduced through an extra force 
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Fig. S3. Plot of the potential energy per particle (normalised by e) of a system 
of N dumbbells, subject to an attractive interaction e. Red and blue curves refer to 
active and passive dumbbells respectively Solid line correspond to a passive gel-like 
initial condition as in the text, whereas dashed lines correspond to a passive gas-like 
initial configuration. 



(F a ) acting on each sphere, with constant magnitude and di- 
rected along the vector joining the front bead of a dumbbell 
to its rear one (front and rear beads are randomly chosen at 
the start of the simulation). We take F r = ^JknT^Rit) with 
R(t) a stationary Gaussian noise with zero mean and variance 
{R(t)R(t')) = 8(t — t'), so that the fluctuation-dissipation the- 
orem holds in the passive limit (F a = 0). Note that Brownian 
dynamics neglects hydrodynamic interactions. 



In our simulations, we set 7 = 2r 



(ma/e) 



1/2 



the time unit, which is set to 1 in our simulations). The to- 
tal simulation time is typically 10 3 r (we chose a time step 
St = 10" 3 r). The value of 7 is relatively large to ensure an 
effectively overdamped motion on the length scale of the par- 
ticle size, which is the relevant regime for the bacteria in our 
experiments. Finally, the magnitude of the active force, F a 
is chosen so as to correspond to a propulsion velocity of ~ 9 
jLtm/s, close to the peak in the velocity distribution in Fig. 3B 
in the main text. The relevant dimensionless number associ- 
ated with the propulsion velocity is FaCr/ksT, which is equal 
to 20 in our simulations and may be thought of as an 'active 
Peclet number' Another dimensionless number relevant to 
our corresponding state theory is the analogous of the pre- 
viously defined quantity /, which can be computed here as 
F a (r c — r m in) /ksT ~ 1.6, where r m in = 2 1//6 cr is the point at 
which the Lennard- Jones potential is minimum. 

Stability of active and passive gels and dependence on the 
initial condition. In order to study aggregation of the dumb- 
bells into rotating clusters (cf. Fig. 2B in the main text), we 
started the simulations from an initial configuration in which 
particles are randomly positioned at a (number) density of 
p = 5 x 10 -3 cr -3 , corresponding to a concentration of about 
5 x 10 9 cells/ml. To compare with confocal micrographs, Fig. 
2B shows a 2cr-thick slab of the simulation box. 

In order to study the effect of activity on phase separa- 
tion (cf. Fig. 4 in the main text), we prepared an initial 
configuration with a gel-like structure (p = O.lcr -3 ). We then 
performed two series of simulations, considering first a system 
of passive dumbbells, and then a solution of self-propelled ones 
(with F a chosen as detailed in the previous section). All the 



3 This is because F a cr / k B T — vcr/D where v is the propulsion speed and D — k B T/~f is 
the diffusion coefficient of a passive sphere of size a. 



Footline Author 



PNAS I Issue Date | Volume | Issue Number | 9 



simulations in both series started from the same gel-like ini- 
tial condition, and in both cases we varied e, which controls 
the depth of the depletion attractive potential, which in ex- 
periments depends on the concentration of NaPSS. We then 
measured in each run the potential energy per particle so as 
to estimate the minimum value of e which is required to not 
melt the initial gel. Figure 4 in the main text shows that such 
minimal values are equal to ~ 22 for the passive case and to 
~ 68 for the active one. 

To estimate the effect of the initial configuration, we per- 
formed a second batch of BD simulations, again considering 
both passive and active systems separately, but now starting 
from a gas of passive dumbbells. The resulting curves for the 
potential energy per particle as a function of e relative to the 
two initial conditions are shown together in Fig. S3. In the 
passive case, we find no effect of the starting configuration; 
but in the active case we find hysteresis. When the initial 
configuration is a gas, a stronger attraction is needed to form 
a gel than to keep a gel stable. 

Finally, we have repeated these simulations (starting from 
a passive gel) for a range of densities (0.005a -3 < p < 
0.25cr -3 ) and computed the critical values of e/fcsT, which 
we call x*, after which the gel is stable. In all cases, x* is 
larger in the active case. The ratio between the active and 
passive x*'s depends very weakly on the density of the system 
throughout the range we have simulated (cf. the theoretical 
predictions shown in Fig. S2), and closely approaches ~ 3, as 
also found in our experiments. 



Torque exerted by a bacterium on the cluster 

Here we estimate the torque exerted by a single bacterium on 
the surface of a cluster, Eq. (4) of the main text. 

For convenience, we assume that the flagellum is oriented 
tangent ially to the cluster. We will further neglect that the 
flagellum is a helix, since its radius is typically much smaller 
than its length (~ 0.2/mi vs ~ 10/xm [9]), and replace it by a 
thin cylinder of length I. Note that helicity is needed to gen- 
erate a non-zero -Z^prop 

but results only in a small correction 

to what follows. 

The force generated by a small portion of the flagellum of 
length dz can be approximated by F prop dz/l, while a force of 
the same magnitude and opposite direction is applied locally 
to the fluid. The former force propagates along the flagellum 
and applies a torque on the cluster of magnitude F prop Rdz/l, 
where R is the radius of the cluster. The latter force, applied 
by the flagellum to the fluid, generates a hydrodynamic flow 
that results in a torque of the opposite direction applied to 
the cluster. For a force applied to the fluid at some distance z 
along the flagellum, see Fig. S4, the magnitude of this torque 
is given by [12] 



dT hyd = 



R 



[21] 



where A is the distance along the line connecting the centre 
of the cluster and the point where the force is applied, and f± 
is the component of the applied force F prop dz/l perpendicular 



to that line. From simple geometry we have 
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The total torque applied to the cluster is thus the sum of the 
"direct" and "hydrodynamic" torques 



R 



1- t= = , [231 



which is Eq. (5) of the main text. 



Cluster sizes 

Note that a significant range of cluster sizes is represented in 
the data presented in Figure 5 of the main text: If we had 
used the alternative (equally valid) variable of the number of 
bacteria in a cluster, N — V/v — (i?/cr) 3 , where V and v are 
the volumes of the cluster and an individual bacterium, and 
R and a ~ 1/xm are their corresponding radii, then our x-axis 
would have spanned 2 to ~ 50. 

Two factors probably control the size of the biggest clus- 
ters we observe. First, in the equivalent equilibrium system, 
the size distribution of pre-transition clusters is exponential 
(see refs. [9,10] in the main text), with the average increas- 
ing as one gets closer to the phase boundary. The success of 
our quasi-equilibrium approach suggest that this may also be 
true in our active system. On the other hand, the flagella of 
the surface bacteria create a corona around the cluster that 
prevents other bacteria from joining it. At small cluster sizes, 
this corona will be very sparse, while at larger R most of the 
space around the cluster will be occupied by the flagella thus 
screening the cluster from other bacteria and preventing it 
from growing beyond some critical size. 




Fig. S4. Sketch of the geometry used to compute the torque generated by a bac- 
terium on the surface of a cluster to which it belongs. The cluster is represented as 
a sphere of radius R. We also show the force dipole (dz/) applied at the cell body 
tangentially to the cluster and at some distance z along the flagellum. 
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